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The polaron binding energy Ep in undoped parent cuprates has been determined to be about 1.0 eV 
from the unconventional oxygen-isotope effect on the antiferromagnetic ordering temperature. The 
deduced value of Ep is in quantitative agreement with that estimated from independent optical data 
and that estimated theoretically from the measured dielectric constants. The substantial oxygen- 
isotope effect on the in-plane supercarrier mass observed in optimally doped cuprates suggests that 
polarons are bound into the Cooper pairs. We also identify the phonon modes that are strongly cou- 
pled to conduction electrons from the angle-resolved photoemission spectroscopy, tunneling spectra, 
and optical data. We consistently show that there is a very strong electron-phonon coupling feature 
at a phonon energy of about 20 meV along the antinodal direction and that this couphng becomes 
weaker towards the diagonal direction. We further show that high-temperature superconductivity 
in cuprates is caused by strong electron-phonon coupling, polaronic effect, and significant coupling 
with 2 eV Cu-O charge transfer fluctuation. 



I. INTRODUCTION 

Developing the microscopic theory for high-Tc super- 
conductivity is one of the most challenging problems in 
condensed matter physics. Eighteen years after the dis- 
covery of the high-Tc cuprate superconductors by Bed- 
norz and Miiller^, there have been no microscopic the- 
ories that can describe the physics of high-Tc super- 
conductors completely and unambiguously. Due to the 
high Tc values and the observation of a small oxygen- 
isotope effect on Tc in a 90 K cuprate superconduc- 
tor YBa2Cu307_y (YBCO)^~^, many theorists beHeve 
that the electron-phonon interaction is not important in 
bringing about high-Tc superconductivity. Most physi- 
cists have thus turned their minds towards alternative 
pairing interactions of purely electronic origin. 

On the other hand, there is overwhelming evi- 
dence that electron-phonon coupling is very strong 
in the cuprate superconductors'"^'^. In partic- 
ular, various unconventional oxygen-isotope effects 
Zhao and his coworkers have observed since 1994 
clearly indicate that the electron-phonon interactions 
are so strong that polarons/bipolarons are formed 
in doped cuprates'"^'^^'"'^^'^''"^^'^^ and manganites^^^^^, 
in agreement with a theory of high-temperature 
superconductivity^'^ and the original motivation for 
the discovery of high-temperature superconductivity^. 
However, such clear experimental evidence for strong 
electron-phonon interactions from the unconventional 
isotope effects has been generally ignored. In the 2001 
Nature paper^^, Lanzara et al. appear to provide ev- 
idence for strong coupling between doped holes and 
the 70 meV half-breathing phonon mode from angle- 
resolved photoemission spectroscopy (ARPES). They 
further show that this 70 meV phonon mode can lead 
to d-wave pairing symmetry and is mainly responsible 
for high-temperature superconductivity^"*. Very recently, 
Devereaux et al.'^^ have proposed that the 40 meV Big 



phonon mode rather than the 70 meV half-breathing 
phonon mode is responsible for d-wave high-temperature 
superconductivity. This pairing mechanism contradicts 
the very recent ARPES data, which show that multiple 
phonon modes at 27 meV, 45 meV, 61 meV, and 75 meV 
are strongly coupled to doped holes in deeply underdoped 
La2_a;Sra;Cu04^''. The strong coupling to the multiple 
phonon modes is not in favor of d-wave gap symmetry 
but may support a general s-wave gap symmetry 

Here we determine the polaron binding energy Ep 
for undoped parent cuprates from the unconventional 
oxygen-isotope effect on the antiferromagnetic ordering 
temperature (Tat)'. The determined value (^1.0 eV) 
of Ep is in quantitative agreement with that estimated 
from independent optical data^^ and that estimated the- 
oretically from the measured dielectric constants^^. The 
substantial oxygen-isotope effect on the in-plane super- 
carrier mass observed in optimally doped cuprates^'*^'*^ 
indicates that polarons are bound into the Cooper pairs. 
We also identify the phonon modes that are strongly cou- 
pled to conduction electrons from the angle-resolved pho- 
toemission spectroscopy, tunneling spectra, and optical 
data. We consistently show that there is a very strong 
electron-phonon coupling feature at a phonon energy of 
about 20 meV along the antinodal direction and that this 
coupling becomes weaker towards the diagonal direction. 
We further show that high-temperature superconductiv- 
ity in cuprates is caused by strong electron-phonon cou- 
pling, polaronic effect, and significant coupling with 2 eV 
Cu-O charge transfer fluctuation. 

II. OXYGEN-ISOTOPE EFFECT ON Tn IN 
LA2CUO4 

The antiferromagnetic order (AF) observed in the par- 
ent insulating compounds like La2Cu04 signals a strong 
electron-electron Coulomb correlation. On the other 
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hand, if there is a very strong electron-phonon coupling 
such that the Migdal adiabatic approximation breaks 

down, one might expect that the antifcrromagnetic ex- 
change energy should depend on the isotope mass. Fol- 
lowing this simple argument, Zhao and his co-workers 
initiated studies of the oxygen isotope effect on the AF or- 
dering temperature in several parent compounds in 1992. 
A noticeable oxygen-isotope shift of T/v was consistently 
observed in La2Cu04^. 
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FIG. 1. The temperature dependence of the susceptibility 
for the and ^^0 samples of the undoped La2Cu04 (upper 
panel), and of the oxygen-doped La2Cu04-Hj/ (lower panel). 
After^ 

Fig. 1 shows the temperature dependence of the 
susceptibility for the ^^O and ^^O samples of un- 
doped La2Cu04 (upper panel), and of oxygen doped 
La2Cu04-|-j, (lower panel). One can see that the AF or- 
dering temperature Tjq for the sample is lower than 
the ^^O sample by about 1.9 K in the case of the un- 
doped samples. For the oxygen-doped samples, there is 
a negligible isotope effect. 

It is known that the antifcrromagnetic properties of 
La2Cu04-i-3/ can be well understood within mean- field 
theory which leads to a T/v formula^°: 



(1) 



where J' is the intcrlayer coupling energy, ^{Tj^) is the 
in-plane AF correlation length at Tn with £,{Tn) oc 
exp(J/T;v) for y = (J is the in-plane exchange en- 
ergy). When Tj\f is reduced to about 250 K by oxygen 



doping, a mesoscopic phase separation has taken place so 
that ^{Tn) = L (Ref.^^), where L is the size of the an- 
tifcrromagnctically correlated clusters, and depends only 
on the extra oxygen content y. In this case, we have 
Tn = J'{L/a)^. Since L is independent of the isotope 
mass, a negligible isotope shift of T/v in the oxygen-doped 
La2Cu04-i-y suggests that J' is independent of the isotope 
mass. Then we easily find for undoped compounds 
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where B = 2J/Tn ~ 10. From the measured isotope 
shift of Tn for the undoped samples, we obtain l^J / J 
~ -0.6%. 

Recently, Eremin et al.^^ have considered strong 
electron-phonon coupling within a three-band Hubbard 
model. They showed that the antifcrromagnetic ex- 
change energy J depends on the polaron binding energy 
Ep due to oxygen vibrations, on the polaron binding 
energy E^"" due to copper vibrations, and on their re- 
spective vibration frequencies wq and cocu- At low tem- 
peratures, J is given by^^ 



J= Jo(l + 
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Here Apd is the charge-transfer gap, which is measured 

to be about 1.5 cV in undoped cuprates. The oxygen- 
isotope effect on J can be readily deduced from Eq. 3: 
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Substituting the unbiased parameters Twjq = 0.075 eV, 
A J/ J ~ -0.6%, = 1.5 eV, and Aujq/ujo = 6.0% 
into Eq. 4, we find that E^ =1.0 cV. The total polaron 
binding energy should be larger than 1.0 eV since E^'^ 
should not be zero. The parameter- free estimate of the 
polaron binding energy due to the long-range Frohlich- 
type electron-phonon interaction has been made for many 
oxides including cuprates and manganites^^ . The total 
polaron binding energy for La2Cu04 was estimated to 
be about 1 eV (Ref.^^), in excellent agreement with the 
value deduced above from the isotope effect. The polaron 
binding energy can be also estimated from optical data 
where the energy of the mid-infrared peak Em in the 
optical conductivity is equal to 2'^Ep (Ref.^^), where 7 
is 0.2—0.3 (Ref.^^). The peak position E^ was found to 
be about 0.6 eV for Lai.98Sro.o2Cu04 (Ref.^^), implying 
that Ep = 1.0—1.5 eV. This is in quantitative agreement 
with the value estimated from the isotope effect. These 
results thus consistently suggest that the polaron binding 
energy of undoped La2Cu04 is about 1 eV. Doping will 
reduce the value of Ep and thus Em due to screening of 
charged carriers. The optical conductivity data indeed 
show that Em = 0.44 eV and 0.12 eV for x = 0.06 and 
0.15, respectively^*. 
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One may argue that the mid-infrared peak could arise 
from magnetic excitations. A sharp peak feature at about 
0.35 cV was seen in the optical conductivity of the un- 
doped YBa2Cu306 (Ref.^'^). This feature was also seen 
in other undoped cuprates^^'^^ , and can be well explained 
by the phonon assisted two magnon excitation'^^. How- 
ever, this sharp feature is very different from a broad 
peak at about 0.6 eV in lightly doped Lai.g8Sro.o2Cu04. 
In particular, the maximum conductivity for the two 
magnon peak is about 1 (f2~^cm~^)"^^^'^^, which is over 
two orders of magnitude lower than that for the broad 
peak in Lai.98Sro.o2Cu04 (Ref.^^). Therefore, the broad 
peak in doped Lai.98Sro.o2Cu04 cannot have the same 
origin as the sharp peak in the undoped system. 

III. OXYGEN-ISOTOPE EFFECT ON THE 
IN-PLANE SUPERCARRIER MASS 

One of the most remarkable oxygen-isotope effects we 
have observed is the oxygen-isotope effect on the pene- 
tration depth^"^'-'^^"^^'^^. We made the first observation 
of this effect in optimally doped YBa2Cu306.93 in 1994^. 
By precisely measuring the diamagnetic signals for the 
-'^^O and -'^^O samples, wc were able to deduce the oxygen- 
isotope effects on the penetration depth A(0) and on the 
supercarrier density n^. It turns out that An^ ~ 0, and 
AA(0)/A(0) = 3.2 % (Ref.S). These isotope effects thus 
suggest that the effective supercarrier mass depends on 
the oxygen-isotope mass. 
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FIG. 2. The doping dependence of the exponent (/?o) 
of the oxygcn-isotopc effect on the in-plane supercarrier 
mass in La2-a;Sra:Cu04. The exponent is defined as 
Po = -dlnmll/dlnMo. The data are from Ref.^'*'^®'^^ 

In fact, for highly anisotropic materials, the observed 
isotope effect on the angle-averaged A(0) is the same 
as the isotope effect on the in-plane penetration depth 
Aab(O). Prom the magnetic data for YBa2Cu306.93, 



Lai.85Sro.i5Cu04, and Bii.6Pbo.4Sr2Ca2Cu30io+j,, we 
found that AAa6(0)/Aa6(0) = 3.2±0.7% for the three op- 
timally doped cuprates^^. Several independent experi- 
ments have consistently shown that the carrier densities 
of the two isotope samples are the same within 0.0004 
per imit cell®'*''^'''. Therefore, we can safely conclude that 
the observed oxygen-isotope effect on the in-plane pene- 
tration depth is caused only by the isotope dependence 
of the in-plane supercarrier mass m**^. Recently, direct 
measurements of the in-plane penetration depth by low 
energy muon-spin-relaxation (LE/xSR) technique^^ have 
confirmed our earlier isotope-effect results. It was found 
that^^ AAab(0)/Aab(0) = 2.8±1.0%. It is remarkable that 
the isotope effect obtained from the most advanced tech- 
nology (LE/iSR)^^ is the same as that deduced from sim- 
ple magnetic measurements^' 

Fig. 2 shows the doping dependence of the exponent 
(/3o) of the oxygen-isotope effect on the in-plane super- 
carrier mass in La2-xSra;Cu04. Here the exponent is 
defined as /3o = — din m*^ /din Mo (where Mq is the 
oxygen mass). It is apparent that the exponent increases 
with decreasing doping, in agreement with the fact that 
doping reduces electron-phonon coupling due to screen- 
ing. The large oxygen-isotope effect on the in-plane su- 
percarrier mass cannot be explained within the conven- 
tional phonon-mediated pairing mechanism where the ef- 
fective mass of supercarriers is independent of the isotope 
mass"^^. In particular, the substantial oxygen- isotope ef- 
fect on m*J; in optimally doped cuprates indicates that 
the polaronic effect is not vanished in the optimal dop- 
ing regime where the BCS-like superconducting transi- 
tion occurs. This suggests that polaronic carriers may 
be bound into the Cooper pairs in optimally doped and 
overdoped cuprates, in agreement with theory'^®'^^''^^. 



IV. STRONG ELECTRON-PHONON COUPLING 
FEATURES ALONG THE DIAGONAL 
DIRECTION 

In conventional superconductors, strong electron- 
phonon coupling features can be identified from single- 
particle tunneling spectra. For high- Tc cuprates, high- 
quality tunneling spectra are difficult to obtain because 
of a short coherence length. Moreover, due to a strong 
gap anisotropy, the energies of the strong coupling fea- 
tures will depend on the tunneling directions. Only if 
one can make a directional tunneling, one may be able 
to accurately identify the electron-phonon coupling fea- 
tures from the tunneling spectrum. On the other hand, 
the observation of the electron self-energy renormaliza- 
tion effect in the form of a "kink" in the band disper- 
sion may reveal coupling of electrons with phonon modes. 
The "kink" feature at an energy of about 65 meV has 
been seen in the band dispersion of various cuprate su- 
perconductors along the diagonal ("nodal") direction^^. 
Prom the measured dispersion, one can extract the real 
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part of the electron self-energy that contains information 
about coupling of electrons with collective boson modes. 
The remarkable progress in the ARPES experiments is 
that the fine electron-phonon coupling structures have 
been revealed in the high-resolution ARPES data of a 
Be surface^". Very recently, such fine coupling struc- 
tures have been also seen in the raw data of the electron 
self-energy of deeply undcrdoped La2-a;Sra;Cu04 along 
the diagonal direction^". Using the maximum entropy 
method (MEM) procedure, they are able to extract the 
electron-phonon spectral density a^F(w) that contains 
coupling features at 27 meV, 45 meV, 61 meV and 75 
meV. The energies of these coupling features are one-to- 
one correspondences to the phonon energies measured by 
inelastic neutron scattering^". These beautiful ARPES 
data and exclusive data analysis^^ clearly indicate that 
the phonons rather than the magnetic collective mode 
are responsible for the electron self-energy effect. 
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FIG. 3. The real part of the electron self-energy along 
the diagonal direction for a slightly overdoped BSCCO with 
Tc = 91 K (0D91K)*i and for La2-a,SrxCu04 with x = 
0.03 (Ref.^°). The energy scale for BSCCO is shifted down 
by Au = 8.0 moV. The solid vertical lines at 27 mcV, 45 
meV, and 61 meV mark the energies of the pronounced 
phonon peaks in the electron-phonon spectral density a'^F{uj) 
of Lai.97Sro.o3Cu04, which is determined from the MEM 
procedure^". The dashed vertical line indicates the energy 
of an extra pronounced phonon peak (75 mcV) in the super- 
conducting LSCO with X = 0.07. It is interesting that the 
coupling feature at 27 meV in LSCO appears to shift to a 
lower energy of about 20 meV in BSCCO. 

Here we will demonstrate that the fine coupling struc- 
tures also appear in the earlier high-resolution ARPES 
data of a shghtly overdoped BSCCO with = 91 K 
(0D91K)^^. Fig. 3 shows the real part of the electron 
self-energy along the diagonal direction for the 0D91K 
sample at 70 K. In the same figure, we also plot the 
real part of the electron self-energy along the diago- 



nal direction for the nonsuperconducting La2-a:Sra;Cu04 
(LSCO) with a;=0.03. We can clearly see the fine struc- 
tures in the raw data of both LSCO and BSCCO. The 
solid vertical lines mark the energies of the pronounced 
phonon peaks (27 meV, 45 meV, 61 meV) in the electron- 
phonon spectral density a'^F{uj) of the nonsupercon- 
ducting Lai.97Sro.o3Cu04, which is determined from the 
MEM procedure^*^. The dashed vertical line indicates the 
energy of an extra pronounced phonon peak (75 meV) 
in the superconducting LSCO with x = 0.07. In order 
for three pronounced fine structures (45 meV, 61 meV, 
and 75 meV) in the self-energy of BSCCO to be aligned 
with those for LSCO, the energy scale for BSCCO has to 
be shifted down by An = 8.0 mcV. This suggests that 
the superconducting gap along the diagonal direction is 
about 8 meV at 70 K. Using the BCS temperature de- 
pendence of the gap, one finds An = 10 meV at zero 
temperature. 

The finite superconducting gap of about 10 meV along 
the diagonal direction is consistent with a general s-wave 
gap symmetry (s -|- g wave) with eight line nodes (g > 
s)^^. This gap symmetry has double gap features at Am 
and A£) in the superconducting density of states^^, which 
can be seen in the single-particle tunneling spectra and 
in the Andrcev reflection spectra. Various break-junction 
spectra suggest that An =9.5 meV and Am = 26 meV 
in an overdoped BSCCO with = 89 K (Ref.^^), Ar, 
= 12±1 meV and Am = 24±1 meV in an overdoped 
BSCCO with = 86 K (Ref.''^), A^ = 7.5-9.0 meV and 
Am = 15-18 meV in heavily overdoped BSCCOs with Tc 
= 62 K (Ref.^^'^''). The Andreev reflection spectrum also 
indicates that Ajj = 13 meV in an overdoped BSCCO 
with Tc = 85 K (Ref.^^), while the other Andreev re- 
flection spectrum shows a gap feature at Am = 25 meV 
(Ref.*^). Further, the Raman data*^ clearly indicate s -|- 
g wave gap symmetry in a heavily overdoped BSCCO 
with Tc = 55 K. The Raman intensities in both Big 
and B2g symmetries increase linearly with energy up to 
1.3Am and can be extrapolated to nearly zero values at 
zero energy. The linear energy dependence of the Raman 
intensity up to 1.3 Am in the Big symmetry could be con- 
sistent with cither clean s -|- g wave supercond uctivity"^^ 
or very dirty d-wave superconductivity with y^rAM(O) 
~ 1.3 Am (where F is the impurity scattering rate)"*^. 
Very dirty d-wave superconductivity would give a signif- 
icant residual intensity at zero energy^^, in disagreement 
with experiment^''. Very dirty d-wave superconductivity 
would also give a T^ dependen ce of the in-plane penetra- 
tion depth below T* = 0.83^/TAm{0) I-IAm ^ 180 
K (Ref.^°), in contradiction with the observed lincar-T 
dependence below 10 K^^. In Ref.^^, the author has con- 
sistently explained all the relevant experiments in terms 
of the s -|- g wave gap symmetry and also shown that 
some phase sensitive experiments apparently supporting 
a d-wave order parameter symmetry do not contradict 
the s -h g wave gap symmetry. 

It is interesting to note that the coupling feature at 
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75 meV is invisible in the deeply underdoped LSCO {x 
= 0.03), but becomes pronounced in the superconducting 
LSCO {x= 0.07)2" BSCCO (0D91K). This is con- 

sistent with the neutron experiments that clearly demon- 
strate that the coupling to the 75 meV half-breathing 
mode increases with increasing doping^". Further, the 
coupling feature at 27 meV in LSCO appears to shift to 
a lower energy of about 20 meV in BSCCO (see Fig. 3). 

V. STRONG ELECTRON-PHONON COUPLING 
FEATURES ALONG THE ANTINODAL 
DIRECTION 

The electron self-energy effect along the antinodal di- 
rection has been studied for several BSCCO crystals 
(0D91K, 0D71K, and OD58K)52 rj.^^ feature in 
the band dispersion or the peak feature in the electron 
self-energy along the antinodal direction is much stronger 
than that along the diagonal direction. This indicates a 
much stronger electron-boson coupling. One of the puz- 
zling issues is that the energies of the boson modes shift 
to a much lower energy (about 20 meV) and are nearly 
independent of doping''^. Fig. 4 shows the boson energy 
as a function of Tc for several ovcrdopcd BSCCO. The bo- 
son energy from the ARPES data is calculated according 
to Eboson = Ekink - Am, where Eunk is the kink energy 
in the band dispersion, which is found to be equal to 
the peak energy in the electron self-energy^^. The above 
relation that is predicted by theory^'^ was also used by 
the authors of Ref.^^ to extract the mode energy. Since 
the antinodal gap Am is found to be very close to the 
peak energy in the energy distribution curve (EDC)^^, 
one can simply take Am being equal to the EDC peak 
energy. It is apparent that the boson energy is about 20 
mcV for heavily ovcrdopcd BSCCOs and about 16 mcV 
for nearly optimally doped BSCCO. The strong coupling 
feature at about 20 meV also agrees with the electron- 
boson spectral density a^Fljjj) deduced from a break- 
junction spectrum of BSCCO (OD93K)^^, as shown in 
Fig. 5. The spectral density clearly shows strong cou- 
pling features at about 20 mcV, 36 meV, 60 meV, and 72 
meV, similar to the features at about 20 meV, 62 mcV, 
and 75 meV along the diagonal direction (see Fig. 3). 
Because the superconducting gap is very anisotropic, the 
excellent match of the phonon energies obtained from the 
tunneling spectrum and ARPES indicates that the tun- 
neling process of this break-junction should be rather di- 
rectional. Comparing the gap size of 23 meV determined 
from the spectrum^''' with the angle dependence of the 
gap^^'2®, we find that the spectrum may mainly probe 
the superconducting density of states along the directions 
of between 10° and 20° from the antinodal direction. 
It is interesting that the strong coupling feature at 36 
meV is only scon in the tunneling spectrum that mainly 
probes the states near the antinodal regime. The 36 meV 
phonon mode should be the oxygen buckling mode (Big) 



that has been shown to couple more strongly to the states 
near the antinodal direction^^. Such a strong angle de- 
pendence of the coupling strength may also occur to the 
other phonon modes such as the 45 meV mode that has 
a stronger coupling along the diagonal direction. What 
is more puzzling is that the band dispersion along the 
antinodal direction has a single kink feature associated 
with the 20 meV boson mode. A very likely explanation 
is that the phonon modes with energies higher than 30 
meV lie below the bottom of the band continuum and 
thus the kink features for these modes disappear^^. Sig- 
nificantly away from the antinodal direction, the energies 
of all the modes are within the band continuum so the 
kink features for these modes will show up. 



30 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

BSCCO (OD) 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 
• ARPES 








o Tunneling - 




1 




- 






< 


1 




1 







55 60 65 70 75 80 85 90 95 
Tc (K) 

FIG. 4. The boson energy as a function of Tc for several 
ovcrdopcd BSCCO. The boson energy extracted from ARPES 
(Ref.'''^) is calculated according to Eboaon = Ekink — Am, 
where Ekink is the kink energy in the band dispersion, which 
is found to be equal to the peak energy in the electron 
self-energy®^. One data point (open circle) is from the tun- 
neling data (see Fig. 5). 

The much stronger coupling to the 20 meV phonon 

modes along the antinodal direction is rather unusual. 
This is possible if the extended van Hove singularity is 
about 20 meV below the Fermi level and the clcctron- 
phonon matrix element for the 20 mcV phonon modes has 
a maximum around q — Q, where (f is the phonon wavevec- 
tor. The large density of states at the van Hove singular- 
ity (20 mcV below the Fermi level) and strong Fermi sur- 
face nesting along the antinodal direction greatly enhance 
the phase space available for 20 mcV small-if phonons to 
scatter quasiparticles from the states near the antinodal 
regime to the extended saddle points. The first princi- 
ple calculation'"'® indeed shows that unusual long-range 
Madelung-likc interactions lead to very large matrix el- 
ements especially for zone center modes {q = 0), which 
arc mainly related to vibrations of cations (e.g.. La, Sr, 
Ba, Ca). The phonon energies for the vibrations of the 
cations are between 15 meV to 25 meV (Ref.^^). 
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Since the strong coupling feature almost disappears 
above Tc, the 20 meV bosonic mode may be related to 
the magnetic resonance mode rather than the phonons. 
First, the energy of the bonsic mode is nearly indepen- 
dent of Tc (see Fig. 4) while the resonance energy is pro- 
portional to Tc (Rcf."''^). This implies that the 20 mcV 
boson is not the magnetic resonance mode. Second, the 
disappearance of the coupling feature above Tc does not 
necessarily mean that the coupling to the phonon modes 
is irrelevant. Very recent calculation^^ shows that the 
feature of coupling to the Big phonon mode (36 meV) 
is very weak in the normal state. As pointed out by 
these authors^^, the dramatic temperature dependence 
arises from a substantial change in the electronic occu- 
pation distribution and the opening of the superconduct- 
ing gap. In the normal state, the phonon self-energy is 
a Fermi function at 100 K centered at the phonon en- 
ergy, which results in a thermal broadening of AAksT or 
38 meV, significantly larger than the phonon energy (20 
meV). At low temperatures and in the superconducting 
state, the phonon self-energy is sharply defined due to 
the step-function-like Fermi function and a singularity in 
the superconducting density of states. 
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FIG. 5. The electron-phonon spectral density a'^F{ui) for 
a slightly overdoped Bi2Sr2CaCu20s+j/ (BSCCO) crystal, 
which was deduced from a break-junction spectrum' 
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Now we hirthcr show that a strong coupling feature 
in optical data, which was previously explained as due 
to a strong coupling between electrons and the mag- 
netic resonance mode^^'^^, is actually consistent with 
a strong electron-phonon coupling at a phonon energy 
of about 20 meV. It is known that the electron-phonon 
spectral density a'^F{oj) can be obtained through inver- 
sion of optical data. Marsiglio et alS'''^ introduced a di- 
mensionless function W{(jj) which is defined as the sec- 
ond derivative of the normal state optical scattering rate 
t~^{llj) = (flp/4nyR(j^^{oj) multiplied by frequency uj. 
Here ftp is the bare plasma frequency and cr{u)) the nor- 
mal state optical conductivity. Specifically, 



27r duj'^ t{oj) 



(5) 



which follows directly from experiment. Marsiglio et al.^^ 
made the very important observation that within the 
phonon range W{lu) ~ a'^F{Lu). 

In the superconducting state, a phonon mode that is 
strongly coupled to electrons will appear at an energy of 
2A{k) + uipfi (where tOph is the phonon energy), that is, 
the energies of the phonon structures shift up by the pair- 
breaking energy 2A{k)^^. Because the 20 meV phonon 
modes are much more strongly coupled to the states near 
the antinodal regime and because there is a large quasi- 
particle density of states at the maximum gap edge, there 
must be a maximum at 2AM+cOph in W{lv). For slightly 
overdoped BSCCO with = 90 K, Am = 26.0(5) meV 
(Ref.^^'^®), so we should expect a maximum in W{uj) to 
be at about 72 meV. This is in quantitative agreement 
with the result shown in Fig. 6. 
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FIG. 6. The optically determined clcctron-boson spectral 
density W{u!) for a slightly overdoped BSCCO crystal with 
Tc = 90 K. After«». 



Recently, Devereaux et al. have calculated the 
electron-phonon interactions for the oxygen buckling 
mode (Big) and the in-plane half-breathing mode^'''. 
They find that the 36 meV Big mode couples strongly to 
electronic states near the antinodal regime. They use an 
electron-phonon matrix element that is suitable only for 
YBa2Cu307_y where a large buckling distortion occurs. 
For other cuprates, the Cu02 plane is flat and the buck- 
ling effect is negligible. Raman data have indeed shown 
that the coupling constant of the Big mode in BSCCO 
is more than one order of magnitude smaller than that 
in YBCO (Ref.63). Even for YBCO, the coupling con- 
stant of this mode was deduced to be about 0.05 from 
the Raman data®"^, in agreement with the earlier first 
principle calculation^'*. Moreover, if this 40 meV phonon 
were strongly coupled to the electronic states near the 
antinodal regime, one would expect a maximum in VF(w) 
to occur at about 92 meV in slightly overdoped BSCCO 
with Tc = 90 K. This is in disagreement with experiment 



6 



(see Fig. 6). 

Previously, the energy of the maximum in W{uj) was 
claimed to be in quantitative agreement with the theo- 
retical prediction based on the strong coupling between 
electrons and the magnetic resonance mode^^'^^. These 
authors^*'^^ argued that the maximum in W{lu) should 
occiu at about Am + Er , where Er is the magnetic reso- 
nance energy. For BSCCO with Tc = 90 K, Am = 26-28 
mcV and Ej. = 43 meV, so one expects a maximum in 
W{lj) to occur at 69-71 meV, in reasonable agreement 
with the experimental result (see Fig. 6). Later on, more 
rigorous theoretical approach^^ shows that the maximum 
in W{lo) should occur at about 2Am + £'r rather than at 
Am + Er- Then we should expect a maximum in W{u>) 
to occur at 96-98 meV, in disagreement with experiment. 

If the maximum in W{u}) would occur at Am + Er, 
the electron self-energy determined from the optical data 
would also have a maximum at Am + Er. For the op- 
timally doped BSCCO (OP96K) with = 96 K, Am 
= 37.5 meV (Ref.^^), and E^ can be estimated to be 45 
meV using the relation between Er and Tc (Ref.''^). Then 
the maximum in the optically determined electron self- 
energy would occur at 82.5 meV. This predicted value is 
significantly lower than the measured one (96 meV)^^. 

On the other hand, we can quantitatively explain the 
optically determined electron self-energy data in terms 
of the s -|- g wave gap symmetry and electron-phonon 
coupling. Because the superconducting density of states 
have two sharp maxima at Am and Ad for the s -|- g 
gap symmetry^^, the optically determined electron self- 
energy should have two peak features, one at 2 Am + Ei 
and another at 2 Ad + E^, where Ei {= 20 meV) and E^ 
are the averaged mode energies along the antinodal and 
diagonal directions, respectively. From Fig. 3, we find 
that E2 = 53 meV, which is a simple average of the two 
pronounced peak energies (45 meV and 61 meV) in the 
electron self-energy. For OP96K, we predict two peak 
features at 95 meV and 73 meV in the optically deter- 
mined electron self-energy. The predicted 73 meV peak 
feature will show up as a shoulder below the dominant 
peak feature at 95 meV. For an overdoped BSCCO with 
Tc = 82 K, Am = 22-25 meV (Ref.^^.es^ ^^^^ ^ 12 
meV (see above), so two peak features will be located at 
64-70 mcV and 77 meV, respectively. For an overdoped 
BSCCO with Tc = 60 K, Am = 14 meV and A^ = 9 meV 
(Ref.^^), so two peak features will show up at 48 meV and 
71 meV, respectively. All these predicted peak features 
are in quantitative agreement with experiments^. 

It is interesting to note that the spectral density shown 
in Fig. 5 is extracted from a break-junction spectrum that 
has a clear dip feature at an energy of about 47 meV 
above the gap"'''"'. The similar dip features are also seen 
in the ARPES spectra along the antinodal direction^^. 
In fact, the dip features of the superconducting density 
of states occur approximately at the valley energies and 
cut-ofi^ energy of the spectral density"^ For exam- 
ple, the cut-off energy of the spectral density for Pb is 



about 9 meV (see Fig. 14 of Ref.^^) and the dip feature 
also occurs at 9 meV (see Fig. 14 of Ref.^^). On the 
other hand, if there is a single sharp peak in the spectral 
density, the dip feature will be slightly above the mode 
energy^^. This implies that the dip energy measured from 
the gap is the upper limit of the mode energy. Because 
the peak features in the spectral density can be broad- 
ened by the strong coupling effect and disorder, the dip 
feature should shift to a higher energy towards the under- 
doping region where the coupling is much stronger®'^. In 
the heavily overdoped region, the coupling is weak^'' and 
the peak width is narrow, so the dip energy is close to 
the mode energy. This can naturally explain why the dip 
energy of the UD70K sample is about 50 meV (Ref.^^), 
while the dip energy of the OD62K sample is about 24 
meV (Ref.^^), which is slightly greater than the mode 
energy deduced from the electron self-energy above. In 
contrast, if one assumes that the dip energy measured 
from the gap is equal to the magnetic resonance energy, 
one cannot selfconsistently explain the dip energy of 50 
meV in the UD70K sample and of 28 meV in the OD87 
K sample^^. The magnetic resonance energies for the 
UD70K and OD87K samples should be about 30 meV 
and 40 meV, respectively^^. 



VI. STRONG COUPLING BETWEEN 
ELECTRONS AND CU-O CHARGE-TRANSFER 
EXCITATION. 
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FIG. 7. The superconducting to normal-state reflectance 
ratio, Rs/Rn, for (BiPb)2Sr2Ca2Cu30io with Tc = 105 K. 
The figure is reproduced from'^^. 

In addition to strong electron-phonon interactions, 

there is a pronounced coupling feature at an energy 
of about 2 eV in the optical reflectance data^"''^'"'. In 
Fig. 7, we plot the superconducting to normal-state re- 
flectance ratio, Rs/Rn for (BiPb)2Sr2Ca2Cu30io. It is 
clear that a strong coupling feature appears at about 
2 eV. A similar strong coupling feature was also seen 
in YBCO, TlaBaaCaaCusOio, and TlsBazCaiCuaOg 
(Ref.''^). Both the temperature and energy dependence 
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of the optical structure can be well described within 
Eliashberg theory with an electron-boson coupling con- 
stant of 0.30-0.35 (Rcf.^^). Because the energy scale of 
this bosonic excitation is similar to the Cu-0 charge 
transfer gap, it is likely that this high-energy bosonic 
mode corresponds to the Cu-0 charge transfer excita- 
tion. A recent calculation based on three-band Hub- 
bard model has indeed shown that the high energy Cu- 
O charge fluctuation can lead to a significant attractive 
interaction between conduction electrons and that the 
pairing symmetry is of extended s-wave {AigY^. The ex- 
tended s-wave pairing symmetry is consistent with the 
conclusion^^ drawn from the comprehensive data analy- 
ses on nearly all the experiments that are used to test 
the gap symmetry. 

VII. PAIRING MECHANISM 

In two of our previous papers^^'^'', we have proposed 
the pairing mechanism for optimally doped and over- 
doped cuprates. The long-range Frohlich-type electron- 
phonon interaction and the short-range interaction of 
electrons with high-energy phonons lead to the formation 
of polarons. These interactions along with the coupling 
of electrons to the high-energy electronic excitations pro- 
duce a negative value for the effective Coulomb pseu- 
dopotential /i* . The polarons arc bound into the Cooper 
pairs due to the negative fi* and additional attractive 
interaction caused by the retarded electron-phonon in- 
teraction with the 20 meV phonon modes. The prob- 
lem could then be solved within Eliashberg equations 
with an effective electron-phonon spectral density for the 
low-energy phonons and a negative Coulomb pseudopo- 
tential produced by the high-energy phonons and other 
high-energy bosonic excitations of purely electronic ori- 
gin. Within this simplified approach, we are able to con- 
sistently explain the observed negligible isotope effect on 
Tc, substantial isotope effect on the supcrcarrier mass, 
large reduced energy gap, and high T,, value. 

In this modified strong-coupling model, the effective 
electron-phonon coTipling constant X^p for the low-energy 
phonons is enhanced by a factor of fp = exp(g^). Here 
= A/ljh, ^ is a constant, and ujh is the frequency 
of the high-energy phonon mode^^. The value of can 
be evaluated from the mid-infrared optical conductivity 
which exhibits a maximum at Em ~ 0.12 cV for optimally 
doped BSCCO^^ and LSCO^*. With Em = 0.12 eV, Hluh 
= 75 meV, we find g'^ = Em/ {2huiH) = 0.8, leading to fp 
= 2.2. 

From the spectral density shown in Fig. 5, we can ex- 
tract the effective electron-phonon coupling constant Xep 
for the low-energy phonon mode, that is, X^p ~ 2.6. If 
there were no polaronic mass enhancement, the coupling 
constant contributed from the low-energy phonons would 
be 2.6//p= 1.2. With n* = 0.1 and Xep = 1.2, we calcu- 
late Tc = 18 K according to a Tc formula^^ 



ksTc = 0.25hV< w2 >[exp(2/Ae//) - l]"^/^ (6) 
where 

Xeff = {Xep - M*)/[l + + Xep/J-*t{Xep)], (7) 

The function t{Xep) is plotted in Fig. 2 of Ref.^^. In the 

present case, ixP- > is contributed only from the low- 
energy phonons and equal to 20 meV. Therefore, without 
the polaronic effect, would not be higher than 20 K. 
On the other hand, with the polaronic effect, X^p = 2.6 
and 11* may be close to zero, leading to a Tc of about 54 K. 
Thus the polaronic effect enhances Tc significantly, but 
electron-phonon coupling alone cannot explain supercon- 
ductivity above 100 K in optimally doped cuprates unless 
the polaronic effect can make n* < -0.15. 

In order to explain superconductivity above 100 K, it 
may be essential to consider the coupling to the high- 
energy electronic excitations. One of the high-energy ex- 
citations is the Cu-0 charge transfer excitation at about 
2 eV, as seen from the optical experiments^^'^^. Since 
the high-energy phonon modes couple to electrons nona- 
diabatically, it is likely that the coupling of electrons to 
the 2 eV boson mode should also be nonadiabatic. Here 
we simply ignore the nonadiabaticity, and consider two 
(5-functions in the electron-boson spectral density to es- 
timate Tc- One (5-function is at huji = 20 meV with 
the coupling constant Ai = 2.6, another at Jiu)2 = 2100 
meV with the coupling constant A2 = 0.3. By solving 
the s-wave Eliashberg equations with the above spectral 
density, we find that Tc = 106 K for fj,* = 0.1. In order 
to obtain Tc = 105 K with one ^- function, we need the 
following parameters: Ai = 2.9, huji = 39.5 meV, and /i* 
= 0.1. If we take /i* = 0, we need Ai = 2.9 and fiuji = 
33.6 meV to have Tc — 105 K. Because we have used the 
(5-functions, the calculated Tc values should be the upper 
limits. 



VIII. CONCLUSION 

We have determined the polaron binding energy Ep 
for undoped parent cuprates from the unconventional 
oxygen-isotope effect on the antiferromagnetic ordering 
temperature (Tjv). The deduced value (about 1.0 eV) 
of Ep is in quantitative agreement with that estimated 
from independent optical data and that estimated theo- 
retically from the measured dielectric constants. The po- 
laron binding energy should be large enough to overcome 
the intersite Coulomb interaction to form intersite bipo- 
larons in deeply underdoped cuprates, in agreement with 
theory and experiment^^'^^. The substantial oxygen- 
isotope effect on the in-plane supercarrier mass observed 
in optimally doped cuprates suggests that polarons are 
bound into the Cooper pairs. The bipolaron picture may 
be irrelevant for optimally doped cuprates because the 
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superconducting transition is of BCS-like. We also iden- 
tify the phonon modes that are strongly coupled to con- 
duction electrons from the anglc-rcsolvcd photocmission 
spectroscopy (ARPES), tunneling spectra, and optical 
data. We consistently show that there is a very strong 
electron-phonon coupling feature at a phonon energy of 
about 20 meV along the antinodal direction and that this 
coupling becomes weaker towards the diagonal direction. 
We further show that high-temperature superconductiv- 
ity in cuprates is caused by strong electron-phonon cou- 
pling, polaronic effect, and significant coupling with 2 eV 
Cu-0 charge transfer fluctuation. The role of the antifer- 
romagnetism in superconductivity appears to be insignifi- 
cant because the antiferromagnctic and superconducting 
phases do not mix, as shown recently by Bozovic and 
coworkers^^. 
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